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Abstract. We place limits on semiclassical fluctuations that might be present in the primor- 
dial perturbation spectrum. These can arise if some signatures of pre-inflationary features 
survive the expansion, or could be created by whatever mechanism ends inflation. We study 
two possible models for such remnant fluctuations, both of which break the isotropy of CMB 
on large scales. We first consider a semiclassical fluctuation in one Fourier mode of primordial 
perturbations. The second scenario we analyze is a semiclassical Gaussian bump somewhere 
in space. These models are tested with the seven-year WMAP data using a Markov Chain 
Monte Carlo Bayesian analysis, and we place limits on these fluctuations. The upper bound 
for the amplitude of a fluctuation in a single Fourier mode is a < 10~^, while for the Gaussian 
bump a < 10~^. 
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1 Introduction 

The standard model of cosmology assumes a homogeneous and isotropic universe, and this as- 
sumption is in good agreement with the observational data from the CMB (cosmic microwave 
background) radiation and galaxy surveys [1]. The standard ACDM model invokes a period 
of exponential expansion, and the observed perturbation spectrum arises from quantum fluc- 
tuations in the free-field theory vacuum (Bunch-Davies vacuum) of the inflaton field that are 
amplified by the expansion of the universe [2]. At present, we do not have a well-defined 
(renormalizable) quantum field theory for inflation which explains not only the period of 
exponential expansion, but also how it ends. It is important to test how well the ACDM 
scenario works, so that one can constrain field theory models for inflation. 

Possible deviations from standard cosmology have been studied in the literature [3-7]. 
It is possible that some anomalies on large scales have evolved directly from the primordial 
perturbations. If inflation does not last long enough to erase the initial transient contributions 
to the background dynamics, some observable features may be left in the CMB [8]. Another 
possibility of modifying the initial power spectrum is by having inflation that does not start 
in the Bunch-Davies vacuum state (see, e.g. [9, 10]). Multi-field inflationary models may also 
generate non-isotropic perturbations [11]. The authors of [12] have discussed the possibility 
of a linearly modulated primordial power spectrum. They found an improvement in of 
about 9 for 3 extra parameters, using the three-year WMAP data. 

In this paper we study the possibility that the primordial density perturbation spec- 
trum, in addition to the standard quantum fluctuations, has a semiclassical contribution. 
Semiclassical fluctuations arise naturally in field theory if the initial state is not the quantum 
vacuum, but contains semiclassical field configurations such as topological defects (vortices, 
monopoles, skyrmions, domain walls, etc.) which can be produced in cosmological phase 
transitions. We analyze two different models for such fluctuations, and put bounds on their 
parameters using the seven-year temperature data from the WMAP satellite. The models 
provide a functional form for the fluctuations with a few parameters, which can be constrained 
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by the data. Firstly, we consider a semiclassical fluctuation in one of the Fourier modes of 
primordial perturbations. This model could arise if one momentum mode of the inflaton 
field was not in its quantum ground state. The second model we consider has a Gaussian 
density bump, and corresponds to initial conditions where the scalar field deviates from its 
vacuum value in some region of space. Clearly, we are only sensitive to fluctuations on scales 
comparable to the size of the observed universe. Much larger wavelength fluctuations would 
appear as a constant shift in the mean density of our universe. 

This paper is organized as follows. In Section 2 we compute the CMB fluctuations, 
including additional semiclassical contributions. This result is then implemented in the like- 
lihood calculation for CMB data in Section 3 for the specific fluctuations we study. The 
details of the data analysis are given in Section 4. We present our results in Section 5, and 
we conclude in Section 6. 

2 CMB Perturbations with a Semiclassical Contribution 

In this section, we describe how to compute the CMB fluctuations including an additional 
semiclassical contribution. We first discuss the general case of an arbitrary semiclassical 
contribution, and then restrict ourselves to the two cases of a single Fourier mode, or a 
Gaussian bump, which are analyzed in this paper. 

We start with a quick review of the standard CMB analysis, and then discuss its modifi- 
cation. The initial perturbations can be expressed in terms of the initial gauge invariant cur- 
vature perturbations C(k) on uniform density hypersurfaces. For homogeneous and isotropic 
Gaussian perturbations the initial curvature perturbations are completely described by the 
two-point correlation function, the power spectrum 

(C(k)C*(k')> = (2vr)3<53(k - k')P(fe) . (2.1) 

Higher order correlations of C can be given in terms of the two-point correlations if the 
fluctuations are Gaussian. The power spectrum P{k) is the usual scale invariant Harrison- 
Zel'dovich spectrum if the inflationary theory is a free scalar field. 

The CMB observable is the temperature anisotropy 0(x, n), the temperature fluctuation 
in direction n as seen by on observer located at x.^ It is more convenient to Fourier transform 
the X dependence, and use 0(k, n). The temperature anisotropies 0(k, n) can be expressed 
in terms of the initial gauge invariant curvature perturbations C(k) 

e(k,n) = C(k)^^^^ (2.2) 

where the ratio Q{k,h- h)/(^{k) does not depend on the initial curvature perturbations. It is 
determined from the evolution of G and and only depends on the magnitude of k and the 
direction of fi relative to k [13]. The angular dependence of the temperature anisotropies is 
written in terms of the spherical harmonics, 

e(x, n) = ^ ai^{x)Yi^{h) (2.3) 

Im 

and 

aU^) = I 0^^'"'' I ^:n(A)0(k, n) ■ (2.4) 

^We only have data when x = xo, our location in the universe. We choose coordinates so that xo = 0. 
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The observed CMB map 0(xo,n) allows one to measure the two-point temperature 
correlation 

Mimi'm' = (az„(0)<,„,(0)) . (2.5) 

All higher order temperature correlations are given in terms of this if the fluctuations are 
Gaussian. 

Expanding 0(A;, k • n) into Legendre polynomials 

@{k, k • n) = Y,{-iy{2l + l)PiCk ■ n)@i{k) (2.6) 

the covariance matrix (2.5) takes the form 

Mlml'm' = ^ll'^mm'Cl (2.7) 

with 



Ci = ^Jdk k^P{k)\Ri{k)f 



H,W = ff. (2.8) 

which is the theoretical prediction for the CMB power spectrum, and can be compared with 
observations. 

In our models, the fluctuation can be written as 

C(x) = Cc(x) + Cci(x) (2.9) 

where Q and (^i are the quantum and classical components of C. Such a form arises naturally 
in a quantum field theory; for example, the quantum field around a semiclassical background, 
such as a soliton, has exactly the same decomposition as in Eq. (2.9), (/'(x) = (/>g(x) + (/)ci(x), 
where (/>ci(x) is the soliton field configuration [14]. 

The usual assumption for the quantum fluctuations is that the one-point average van- 
ishes, 

(C,(x)) = (2.10) 

and the two-point correlation is as in Eq. (2.1). The additional semiclassical component Cci 
is an ordinary function, not a quantum operator, so that 

(Ccl(x)) =Ccl(x), 

(Cci(x)C,(y)) = 0, 

(Cci(x)Cci(y)> = Cci(x)Cci(y) . (2.11) 

In Fourier space, this gives 

(C(k)) = Cci(k) , 

(C(k)C(k')> = Cci(k)C*i(k') + (2^)3<53(k - k')P(fc) . (2.12) 
for the initial C spectrum. 
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Combining Eq. (2.9) with the previous equations, one finds that 

4m = M^)) = 1^/ ^^^^^ k^YtMRmU^) (2.13) 
and the two-point correlation is 

Ml^Vm' = {aUmmm) = 4yfm' + ^^ImVm' (2-14) 

where the first term is the new contribution, and the second term is the standard result 
Eq. (2.7,2.8), now denoted by a superscript q to distinguish it from the total M. In general, 

the n-point correlations of the shifted quantities aim — OLlra ^^^^ same value as in the 
theory without a semiclassical fluctuation. This is true even if the correlations are not 
Gaussian. 

To analyze the experimental data, it is also useful to have formulae for the temperature 
correlations in pixel space. The temperature Aj in a pixel in direction rij is 

Ai = j dne(n)Si(n) (2.15) 

where Bi{h) is the beam pattern of pixel i (including the pixel window function). The 
beam pattern is specific to the experiment. Usually the beam patterns have the same shape 
for every pixel and are axially symmetric around the center of the pixel, as is the case for 
WMAP. If we denote the direction to the center of the pixel by fij then the beam pattern 
can be decomposed into spherical harmonics 

B,ih) = Y,BiYUn,)YCm{n), (2.16) 

Im 

and Bi does not depend on the pixel. Then 

Af) ^ {A,) = Y,{2l + l){-iy I ^BiPi{k-n,)Ri{k)Uk), (2.17) 

and the two-point pixel-pixel temperature correlation is 

(A,A,) = Af )Af + (A,A,)^ (2.18) 
where (AjA^)^ is the standard (quantum) contribution 

(A.A,)^ = Ml,^,BiB!,YUni)Yirm'{n,) 

Iml'm' 

= Y,Ci\Bi\^Yimim)Yanj). (2.19) 

Im, 

Equations (2.17,2.18,2.19) form the basis for the likelihood analysis of the next section. 
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3 Likelihood Calculation 



We use the likelihood method to compare the theory with the experimental data. The 
likelihood is calculated in real space. For Np pixels the likelihood function has the form 

where A is the vector of pixels of measured temperature anisotropics and Cij is the covari- 
ance matrix, including noise. The theoretical covariance matrix is obtained by transforming 
Mimi'm' into real space, and is given by Eq. (2.19), 

a, = {A,Aj)^, (3.2) 

to which the noise matrix must be added. The likelihood is a function of the cosmological 
parameters through the dependence of C on the cosmology. 

In the modified theory we consider, there is a one-point function, Eq. (2.17), and a 
modfied two-point function Eq. (2.18). Thus to analyze the data, we can use the modified 
likelihood function 



^= (2.)«./^(detC)V. -n-2[^-^n ^-'[A-A^JJ. (3.3) 

replacing the observed temperature values by their difference from the expected mean, Aj 
Aj — A-'^^^ The covariance matrix C remains the quantum covariance matrix Eq. (3.2) plus 
instrumental noise. Equation (3.3) can be used to analyze any semiclassical fluctuation in 
the primordial universe. In this paper, we restrict ourselves to two interesting functional 
forms, the single Fourier mode, and the Gaussian bump. 

3.1 Single Fourier Mode 

Consider the case of a semiclassical perturbation in a single momentum mode. We choose 

Cci(x) = ao cos(ko -x + a) (3.4) 

which is given in terms of 5 parameters — the amplitude ao, the wavelength A = 27r/k and 
the direction (the 3 parameters in kg), and the phase a at Earth. From Eq. (2.13), 



27rao(-i)' Riiko) YiUk) (e'" + (-l)'e-^") (3.5) 



and 



^(ci) ^ ^^(2/ + l)(-iyBiPi{h, ■ k) (e'" + (-l)'e-^°) Ri{ko) (3.6) 
I 

3.2 Gaussian Fluctuation 

The second case studied, the Gaussian fluctuation, has 

C^^'Hx) = aoe-(''-^^)'/'"' (3.7) 
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which also has 5 new parameters — the amphtude oq, the radius of the Gaussian w, and the 
position of the center Xc. The new contributions are given by 

«Jf (0) = ^aow\^/^{-iy J edke-^'^''^Ri{k)ji{kT)Y,Uv) (3.8) 

where we have used r = Xc. 

Af) = -l=ao J](-l)'(2/ + l)i?^PKA*-i-) / k^dke-^''^Ri{k/w) ji{kr/w). (3.9) 

We use Eq. (3.3) for the Hkehhood function in our data analysis. It depends on the cos- 
mological parameters through C, and the parameters in the semiclassical fluctuation through 
A^^^). In principle, one should maximize the likelihood with respect to all the cosmological 
parameters as well as the five new parameters. We will use a simpler approach. Since we 
are interested in the breaking of isotropy on large scales, we will assume that the periodic 
fluctuation wavelength is of the same order as the distance to the last scattering surface 
Lq = 14.4 Gpc. This means that our modification will affect only the low-/ part of the co- 
variance matrix. Since the standard cosmological parameters are determined from the whole 
spectrum of /, we fix their values at their best-fit values as given by the seven-year WMAP 
data [15], and only vary the new parameters.^ 

4 Data Analysis 

We use the publicly available CAMB code [16] for calculating the standard covariance matrix 
and Qi{k)/C{k), and a modification of the likelihood code provided by WMAP [15, 17, 18] for 
calculating the likelihood. Since our modification affects only the low-Z part of the spectrum, 
we use the low-resolution part of the likelihood code. The details of the low-/ likelihood 
calculation for WMAP data can be found in [19]. Let us summarize the important points. 
The sky map used for low-/ analysis is the ILC (Internal Linear Combination) map. The 
map is smoothed to 9.183° and degraded to Ngi^e = 16 using the HEALPIX software [20]. 
Smoothing is done to reduce the noise in high /. The map is further masked with the Kp2 
mask, leaving 2482 pixels out of 3072. Foreground contamination, as well as instrumental 
noise in the unmasked region should be negligible [19]. Although it is not recommended to 
use the ILC map for high-resolution analysis, the smoothed and degraded ILC map is reliable 
for low-resolution analysis and is the default sky map used in low-resolution likelihood code 
by WMAP. For comparison, we also do the analysis on the V, W, and Q frequency bands as 
well, smoothed and degraded as described above. 

It is hard to accurately estimate the noise covariance matrix of the smoothed and 
degraded maps since the noises in different pixels are highly correlated. A 1 /iK white noise 
is added to each pixel of the ILC map, and a corresponding term to the diagonal elements of 
the covariance matrix, to aid numerical regularization of the matrix inversion [19, 21]. This 
is small enough to have a negligible effect on the data, but large enough to dominate over 
the instrumental noise present in the data. For the high resolution V, W, and Q maps the 
instrumental noise in each pixel can be estimated by cr = (Tq / V^obs with ctq = 3.137mK (V), 

The values of the cosmological parameters that we use are 100f7b/i^ = 2.227, ^Ich^ = 0.1116, f^A = 0.729, 
Us = 0.966, r = 0.085, A|(0.002 Mpc-^) = 2.42 x lO"''. 
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6.549 mK (W), and 2.197mK (Q) [21]. Noi,s is supplied for each pixel along with temperature 
in WMAP sky maps. In order to estimate the noise in the smoothed and degraded maps 
we generate 1000 simulations of noise at high resolution, then smooth and degrade them. 
The average noise per pixel is estimated in this way to be 0.45 /xK for V, 0.54 fiK for W, and 
0.37 fiK for Q band low resolution maps. Therefore, adding a 1/iK noise to each pixel will 
still dominate over the existing noise, so it is safe to analyze these maps the same way as the 
ILC. 

The residual galactic uncertainty is removed by constructing a foreground template A f 
and marginalizing the likelihood function over it. This is done as follows. An extra parameter 
^ is introduced into the likelihood function (3.1) 

The likelihood is then marginalized over ^ to get 



(27r)^p/2(detC)i/2y AjC-iA/ V H ^/'^"'^/ 



The corresponding formula for Eq. (3.3) is Eq. (4.2) with A — t- A — A^'^^). 

The foreground template used in the WMAP likelihood code is obtained by subtracting 
the ILC map from the V band map. Note that the ILC map is constructed by removing all 
known foreground contaminations, so the foreground marginalization for ILC analysis has a 
very small effect on the results. However, it is essential for analyzing the V, W, and Q bands. 

The monopole and dipole (/ = and 1 terms) are removed from the original maps. 
However, after masking the maps, these terms are re- introduced. Thus the likelihood function 
needs to be marginalized over the monopole and the dipole contributions. This is done by 
introducing large variance I = and 1 terms into the covariance matrix. 

Since we are interested in large scale isotropy breaking, we restrict our analysis to I < 30. 
We check that cutting off at Imax = 30 does not have a significant effect on our results for 
the scales we consider. For the fluctuation in one Fourier mode, the I = 15 terms in Eq. (3.6) 
are already 3 orders of magnitude smaller than the I = 2 terms for the wavelengths A under 
consideration. This means that the effect of the cutoff can be completely neglected. For the 
Gaussian bump, the decay with I is not as fast. In Fig. 1 we plot the dependence of 

m=— Z 

on I for a typical scale w = 0.14Lo = 2Gpc and a typical amplitude of ao = 10~^, keeping 
all the other parameters fixed. The / = 30 term is 2 orders of magnitude smaller than the 
I = 2 term, therefore the cutoff cannot have a significant effect in this case either. 

The CMB photons get lensed by the matter present in the universe on their journey 
from the last scattering surface to the observer. However, the WMAP seven-year temperature 
data by itself is not sensitive to lensing on all scales [22]. Since we are only considering large 
scales {I < 30), the lensing effects become completely negligible. This is why we do not 
include lensing in our analysis. 

We analyze the data using the Bayesian approach. Usually for a large number of param- 
eters the Bayesian approach is computationally easier than the frequentist approach. One of 
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Figure 1. For a Gaussian fluctuation with oq = 10 we show in the units of the CMB 
temperature 2.73 K as a function of I. 

the main computationally efficient methods for Bayesian analysis is the Markov Chain Monte 
Carlo (MCMC) method (see, e.g., [23]), which we will use to obtain confidence intervals for 
our parameters in the Bayesian approach. 

5 Results 

For both of the models considered, the data is analyzed using the likelihood function in the 
space of our new 5 parameters. It is easier to deal with — 21n^ instead of the likelihood ^ 
itself. The analyses is sensitive to likelihood ratios only, therefore the overall normalization 
constant does not matter. Equivalently, any constant can be added to — 21n^ without 
changing the likelihood ratios. To make it easier to compare our hypothesis to the standard 
cosmology we choose to normalize the likelihood function such that 



where is the likelihood for the standard cosmology. 
5.1 Semiclassical Fluctuation in One Fourier Mode 

We construct a Markov Chain using the Metropolis-Hastings algorithm [23] to calculate 
the posterior probability distribution functions for our newly introduced parameters. The 
standard cosmological parameters are kept fixed as before. The prior probability distributions 
are chosen as follows. The distribution for direction is uniform on the whole sphere. We vary 
A in the range [O.SLq, 2.6Lo], so we choose the prior distribution to be uniform on that 
interval and otherwise. The prior distribution for a is uniform on the whole range [0, 27r]. 
We choose a Gaussian distribution with a = 10^^ as the prior distribution for the amplitude. 
This <T is chosen small enough not to allow entry into the nonlinear regime but large enough 
to keep the prior distribution uniform in the region of interest. As we will see later, the 
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Figure 2. Checks of our pipeline witli simulations. We show the posterior distribution p(a:;|Data) 
against ao (left), A (middle), and a (right). In the solid black we show the results from a simulation 
without the presence of semiclassical fluctuations. In dotted blue curve we show the results from a 
simulation with an added semiclassical cosine fluctuation with oq = 10 x 10^'^, A = Lq and a = 0. 
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8.14 


14.20 


Q 
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Table 1. Upper bounds for |ao|/(10 ^) for a semiclassical periodic perturbation for different sky 
maps. 

posterior distribution converges to at scales about two orders of magnitude smaller than 
this (T, meaning that the prior distribution can be thought of as uniform with very good 
approximation. We generate one chain of length 10^ and throw away the first 10^ elements 
for burnin. We compare results from a few chains of the same length and burnin time to 
make sure that they agree well. 

To check this method, we first run the analysis on two simulated maps. The first map 
is a standard universe without any semiclassical fluctuation added. The second simulation 
has a semiclassical cosine fluctuation added with oq = 10 x 10~^, X/Lq = 1.0, a = 0. 
The simulations have the same beam function and noise as the real data. The posterior 
probability distributions for the amplitude p(ao|Data), the wavelength p(A|Data), and the 
phase p(q I Data) are shown in Fig. 2. For simulations without the semiclassical perturbations, 
the amplitude distribution peaks at as expected, while the distributions for A and a are 
essentially flat, without any significant peaks. For simulations with the periodic semiclassical 
perturbations, p(ao|Data) has a clear peak very close to the expected value of 10 x 10~^, 
p(A|Data) peaks sharply near 1.0 x Lq, and p(a|Data) has a sharp peak at a = 0. Those two 
simulations show that the MCMC method is working very well. 

We are now ready to analyze the real sky maps. The posterior distribution functions 
for the amplitude, the wavelength, the phase, and the direction (ILC map only) are shown 
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Figure 3. Results from the WMAP data for a perturbation in a single k mode. We plot the posterior 
distribution p(a; I Data) against oq (upper left), A (upper right), a (lower left), (lower right). In the 
lower right panel we show the ILC map. 

in Fig. 3. There is no detection for non-zero amplitude for all the maps. The non-zero peak 
for the W band correspond to less than 2a. The results from different maps agree reasonably 
well. We also obtain upper bounds on |ao| from p(ao|Data), which are summarized in Table 1. 
The limits obtained from different maps agree reasonably well with each other. 

5.2 Semiclassical Gaussian Fluctuation in Space 

We again restrict our analysis to fluctuations on large scales, which means that the parameter 
w needs to be not much smaller than the distance to the last scattering surface Lq = 14.4 Gpc. 
We also need to make sure that the fluctuation has a significant causal contact with the 
last scattering surface, otherwise it will not have an observable effect on the temperature 
fluctuations. The other issue to keep in mind is that if the center of the bump is very close to 
our position then the corrections to the temperature fluctuations will be nearly constant and 
will be absorbed into the constant background temperature. The same thing is true if the 
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Figure 4. Plot of — 21n^ against r for a Gaussian fluctuation, with all the other parameters fixed 
{w ~ O.SLq, ao ~ 10^'^). The red dashed lines correspond to Lq — w and Lq + w. 
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Figure 5. Check of our pipeline for the analysis of a Gaussian bump. We plot the posterior distri- 
bution p(a;|Data) against oq (left), r (middle), and w (right). In solid black curve we show the results 
from a simulation without the presence of semiclassical fluctuations. In the dotted blue curve we show 
the results from a simulation with a semiclassical Gaussian fluctuation added with oq = 0.5 x 10~'^, 
r = 0.9 Lq and w = 0.1 Lq. 



center is not very close to us but w is very large. As an example, we plot the dependence of 
— 21n^ on r with all the other parameters fixed in Fig. 4. In particular, we choose w = O.SLq 
and ao = 10~^. As we can see, it is peaked near Lq, and the dependence becomes very weak 
for r outside the range [Lq — WjLq + w] (red dashed lines). So we will not consider the values 
of r that are outside that range. This is equivalent to the requirement that the la surface of 
the Gaussian fluctuation must intersect the last scattering surface. 

We use the same MCMC algorithm as for the periodic fluctuation to obtain posterior 
probability distribution functions. For the Gaussian fluctuation the likelihood calculation is 
computationally slower, so instead of constructing one long chain we construct 10 chains of 
length 10^, disregarding the first 10^ elements of each chain for burnin. We again compare 
results from a few analyses to make sure that these lengths are sufficient for convergence. 
The prior distribution for the direction is uniform on the whole sphere. We very w in the 
range [O.OSLq, 0.5Lo] and r in the range [Lq — w,Lq + w] for each value of w. The prior 
distribution for these two parameters is uniform in the area allowed (which is a trapezoid) 
and elsewhere. The prior distribution for the amplitude is Gaussian with a = 10~^. 
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Figure 6. Results from the WMAP data for a Gaussian bump. We plot the posterior distribution 
p(x|Data) against qq (upper left), r (upper right), w (lower left), f (lower right). In the lower right 
panel we show the ILC map. 



This is again small enough to not allow significant probabilities in the nonlinear regime but 
large enough to not have a significant impact on the posterior distribution. As we will see, 
the posterior distribution converges to about 2 orders of magnitude earlier, so the prior 
distribution in the region of interest is essentially uniform. 

As for the periodic fluctuation case, we first check the method with simulated maps 
with and without a Gaussian bump added. Again, the simulated maps have the same beam 
function and noise as the real maps. The results are shown in Fig. 5. For a simulation without 
a Gaussian bump, the distribution for the amplitude p(ao|Data) peaks at as expected, while 
j>(r|Data) and p{w\Data) are essentially flat. For the simulations with a Gaussian bump, the 
parameters are oq = 0.5 x 10~^, r = 0.9 Lq, w = 0.1 Lq. The posterior distributions do indeed 
have sharp peaks around the expected values. As explained before, for a Gaussian bump the 
main effect comes from its intersection with the last scattering surface, so we would expect 
similar results for two bumps of the same size and the same distance from the last scattering 
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1.10 



4.24 
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0.68 



1.27 



5.69 



Q 



0.41 



1.14 



3.64 



Table 2. Upper bounds for |ao|/(10 for a scmiclassical Gaussian perturbation for different sky 
maps. 

surface, but one inside one outside. This explains the smaher peak near r = 1.1 Lq. The 
MCMC method detects the Gaussian bump (or the absence thereof) rehably, and now we 
can turn to the results from the real data. 

The results from the analysis of the WMAP data are shown in Fig 6. The results from 
all of the maps agree very well with each other except for the W band. For the ILC, V, 
and Q maps the peaks of the probability distributions lie near the following values of the 
parameters: oq = 2.5 x 10"^, r = 0.94Lo, w = 0.15 Lq, b = 57°, / = 29°. The amplitude 
probability distribution functions indicate that the amplitude is greater than at about 65% 
confidence level. For the W map the probability distributions peak at the following values 
of the parameters: oq = -4.6 x 10"^, r = 0.94 Lq, w = 0.3 Lq, b = 2.4°, / = 45.0°. The 
probability that the amplitude is less than for the W map is 80%. As we can see, for all of 
the maps the indication of a non-zero amplitude is at a level less than 2a. The upper bounds 
for the magnitude of ag obtained from the probability distribution function are summarized 
in Table 2. 

6 Discussion and Conclusions 

We studied the possibility that the inflationary induced primordial perturbations have a 
semiclassical contribution in addition to the standard quantum fluctuations. We gave a 
general formalism for constraining such a primordial semiclassical fluctuation. Then we 
focused on two specific models of semiclassical perturbations, both of which statistically break 
CMB isotropy on large scales. The MCMC method that we used works very well as seen 
by the results from simulated maps with and without semiclassical fluctuations. Although 
there are some indications of a non-zero amplitude for these fluctuations, they correspond to a 
confidence level less than 2a, so our analysis finds no statistically significant evidence for these 
scenarios. We obtained upper bounds on the amplitude of these semiclassical fiuctuations at 
1, 2, and 3 sigma levels. The bounds for the periodic fluctuation are of the order of 10~^, and 
the bounds for the Gaussian bump are of the order of 10"'^. The reason for weaker bounds 
in the Gaussian case is that its effect on the CMB strongly depends on the position and the 
size of the bump. For a bump too far away, the effect on the CMB is very small, while for 
a bump very close to us, the effect is mostly uniform on the last scattering surface and is 
subtracted away with the monopole term before the data analysis. Even for a bump with a 
significant intersection with the last scattering surface, there is a large contribution into the 
dipole term which is, again, subtracted away before the data analysis. 
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We used the ILC, V, W, and Q maps from the WMAP seven-year release for our analysis. 
The results from the different bands, in particular the upper bounds on the amplitude of the 
semiclassical perturbations, agree reasonably well with each other. 
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